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We study the quantum phase transition properties of a three-dimensional periodic 
array of Josephson junctions with charging energy that includes both the self and 
mutual junction capacitances. We use the phase fluctuation algebra between number 
and phase operators, given by the Euclidean group £'2, and we effectively map the 
problem onto a solvable quantum generalization of the spherical model. We obtain a 
phase diagram as a function of temperature, Josephson coupling and charging energy. 
We also analyze the corresponding fluctuation conductivity and its universal scaling 
form in the vicinity of the zero-temperature quantum critical point. 
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There is significant contemporary interest in quantum critical phenomena. Most studies 
have been carried out in two-dimensionsHI. There are several systems where theoretical results 
have been successfully compared against experiment in artificial networksn and homogeneous 
ultrathin filmsH. There has been some but much less work in the three-dimensional (3- 
D) case, although there is both theoretical and experimental interest in this problem, for 
example, in quantum magnetic systems and high temperature superconductors. There is also 
preliminary progress in fabricating quasi-three dimensional Josephson junction arrays (JJA) 
with ultrasmall junctions, in which quantum fluctuations are essentials. There is also interest 
in the classical limit of bulk high-T^ superconductors where the scaling critical properties 
are dominated by thermal fiuctuationsQ. Closely related is also the physics that governs the 
interplay between local and global superconductivity in granular materials, in which disorder 
may also play an important role. In spite of this interest, however, 3-D quantum-capacitive 
JJA have not been investigated in depth yet. Notably, there appear to be no studies on 
3-D JJA close to the T = quantum-critical (QC) point, where the physics is dominated 
by zero-point quantum fluctuations rather then thermal effects. When the superconducting 
islands can sustain at least one Cooper pair the development of global superconducting 
phase coherence depends on the relative strength of the inter-island Josephson coupling Ej, 
as compared to the charging energy Eq = e^/2C, where C is the junction capacitance. In 
the quantum regime the phase-charge interplay is a direct consequence of the Heisenberg 
uncertainty relations between the island phase (pj and the particle number operator Lj = 
id/d(j)j. In this paper we investigate a general quantum-capacitive model for 3-D JJA on 
the simple cubic lattice. We establish the general phase transition boundary and present 
results for the measurable frequency dependent conductivity in the QC regime. We employ 
a novel non mean-field approach based on the proper quantum phase fluctuation algebra, 
by mapping the 3-D JJA model onto an effectively constrained system - a solvable quantum 
spherical model. 

We start by defining a cubic Josephson junction array with superconducting phases <^j 
at the 3-D lattice sites i. The corresponding effective Euclidean action, in the Matsubara 
"imaginary time" r formulation (0 < r < l/ksT = f3, with T being the temperature) is 
5[(/)]=5c[(/)] + 5j[0], where 

SM = 

M<f>] = 11 f d^J-^J {1 - cos[0.(r) - 0,(r)]} . (1) 

Here Sc [0] defines the electrostatic energy, with Cij being the geometric capacitance matrix 
of the array. This matrix is normally approximated, both theoretically and in experimental 
interpretations as: Cij = {Cg + zCm)Sij — CmJ2d^i,j+d, with the vector d running over 
nearest neighbors, with Cg the self-capacitance and Cm the mutual-capacitance between 
nearest neighbors ( z stands for the coordination number). There are more general forms of 
the full capacitance matrix^, but in our analysis the mutual capacitance approximation is 
sufficient. Finally, Sj[(j)] gives the Josephson energy Ej (with Jij = ii^j for |z — j| = |(i| and 
zero otherwise). 

Most analytical works on quantum JJA have employed different kinds of mean-field-like 
approximationsErt3, which are not fully reliable to treat spatial and temporal quantum phase 




fluctuations. Furthermore, as pointed out recentlyllll, the JJA model (|l]) must encode the 
phase fluctuation algebra given by the Euclidean group E2, that involves the commutation 
relations between particle Lj and phase (ladder) operators Pj = e*'^-': [Li,Pj] = —Pj6ij, 
[Li, Pj] = PjSij and [Pj, Pj] = with the conserved quantity (invariant of the E2 algebra) 



p. p.T = p^. I p^. 



(2) 



Thus, the proper theoretical treatment of a quantum JJA must maintain the constraint (^. 
To proceed we write the partition function Z = J [Hi ^0j] e""^'"^]] for the model (|I]) in terms 
of its path integral represent at ionllj, by introducing the auxiliary complex fields i'iir), which 
replace the original ladder operators Pj. To proceed, we substitute the "rigid" E2 constraint 
given in Eq. (|^), by the weaker spherical closure relation j^^iPiP} = 1, which maintains 
(on average) the original condition of Eq. (H). This substitution allows us to formulate the 
problem in terms of an (exactly) soluble quantum spherical (QS) model (see Ref. . |13D. By 
using the Fadeev-Popov method with the Dirac delta-functional, which facilitates both the 
change of integration variables and the imposition of the spherical constraint we obtain: 



Z = J WV^^iV^l^t 5{y^m''-N\e'^^^ 



\{'D(t)^ 



e~^c[^^\{5[^e^lj,-P:m 



X 5['^miPi-Pf{(P)]. 



(3) 



The convenient way to enforce the spherical constraint is to use the functional analog of the 
(5— function representation 5{x) = J^^{dX/2Tc)e''^^, which introduces the Lagrange multiplier 
A(r) thus adding an additional quadratic term (in the t/^— fields) to the action (|l|). The 
evaluation of the effective action in terms of the ip fields may be organized using the loop 
expansion methodHj. To second order in ipi{^) "we obtain the partition function of the 
quantum-spherical model Z = Zqg: 
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(4) 



where 



5qs[V^, \]=Y.f dTdr'{[{J,, + \5,,) 5{r - r') 
+ r,,(r - r')] ^tir)^A^') - iV(5j,A(5(r - r')} . 



(5) 



Here, Fj^ (r — r') is the two-point phase vertex function related to the phase-phase cumulant 
correlation function Wkj{T — r') by 



Y: f dT"Ti,{T - T")Wkj{T" - t') = 5,:,5(r - r'). 

7- •'0 



(6) 



Explicitly, 
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e(o) 

where Zq is the statistical sum of the "non-interacting" system described by the action i5c[^]. 
Since the values of the phases (pi which differ by 27r are equivalent, the path integral can 
be written in terms of the non-compact phase variables Oj{t), defined on the unrestricted 
interval (—00, +00), and by a set of winding numbers {fij} = 0,±1,±2,..., which are 
integers running from —00 to +00 (and physically reflects the discreteness of the chargeli^), 
so that 0j(r) = %(0) + 2Tiin,jT/l3 + ej{T). 

In the N ^ 00 thermodynamic limit the steepest descents method becomes exact; the 
condition that the integrand in Eq.(^) has a saddle point A(r) = Aq becomes an implicit 
equation for Aq: 

l = ^EG'(k,a;,), (8) 

where G-^k, ui) = [Aq - J(k) + 2Ec + ujj/SEc] with uj^ = 27r£//? {£ = 0, ±1, ±2, . . .) being 
the (Bose) Matsubara frequencies and J(k) the Fourier transform of the Josephson couplings 
Jij, respectively. As mentioned above, we next proceed by assuming that Cij has only the 
nearest-neighbor mutual components. For a 3-D simple cubic lattice we obtain for the 
charging energy 

Ec = ie2|C-l|.. 

- ^.c '^-^r'r'!;:"'" KK)K(K-). (9) 

where Eqc = e^ /{2Cs) is the charging energy for the self-capacitive model and K(x) stands 
for the complete elliptic integral of the first kindll3. Furthermore, 
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where 7 = |(1 + zCJCs)/{CJCs). 

As usual in spherical model calculations the phase boundary (i.e. the location of the 
critical points) is determined by Eq.(|^) from the the upper limit of the eigenvalue spectrum 
max{J(k)} = 3Ej, associated with the onset of the phase transition - in the spherical 
model the Lagrange multiplier Aq "sticks" to that value at criticality (Aq = Aq'''* = 3Ej) 
and stays constant in the whole low temperature phaseEJ. Equivalently, at the critical point 
1/G(k = 0,uji = 0) = 0. Introducing the density of states p{E) = JZ^ [d^k/{27Ty]6{E-J{k)) 
we obtain for the critical line 



?,Ej-E 
pj2Ec{3Ej - E) 



where p{E) = Ps=q{E) and 
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(12) 



with ai = max(— 1, —2 — E/Ej), a2 = min(l, 2 — E/Ej); 6(x) is the unit step function and 
Aoix) = 1. 

The current response to an externally applied electromagnetic field is the conductivity a 
that is experimentally measurable. In the context of two-dimensional JJA there are several 
studies of a e.g. at the superconductor-Mott-insulator using 1/A^ expansion and Monte 
Carlo analysis (see, Ref. 18), the coarse-grained approachEJ and an e-expansionE3. 

The standard Kubo formula relates the conductivity to a two-point current-current cor- 
relation function. Applying an external vector potential A modifies the Josephson coupling 
by introducing a Peierls phase factor according to: Jij —>■ Jij exp{2ei/hc J^ A ■ dl). The 
conductivity is obtained as the second derivative of Zqs given in Eq.( ^) with respect to A. 
After performing the derivatives we obtain (for vanishing magnetic field) the longitudinal 



component of (7{uJu) = (^x 
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where Rq = h/Ae^ = GAbkfi is the quantum unit of resistance. We introduced the 
modified density of states p{E) = Z^[d^k/{27rY]sm'^{h)5{E - J(k)), with p{E) = 
I [ps=o{E) — ps=2{E)], where Ps{E) is given in Eq.([T^) with A2{x) = 2x^ — 1. Evaluating the 
summation over Matsubara frequencies and analytically continuing to real frequencies we 
obtain for the real part a' = cr'smg + "^reg of the complex dynamic conductivity (the imaginary 
part a" can be obtained via the standard dispersion relation) 
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where ?7(x) = Ejp{Ejx) and the parameter 5 = Sx/Ej measures the distance from the 
critical point with 6x = Xq — Aq"^'*, and with Uc = \/32EcSx the threshold frequency for 
particle-hole excitations. 

The real part of the conductivity contains two contributions: first, the Drude weight 
cr'singi^), is singular since it is proportional to 6{uj) and the regular finite-frequency con- 
tribution to the conductivity, a[.^^{uj), which is due to the electromagnetic field induced 
transitions to excited states. The singular part in turn is due to the free charge acceler- 
ation. This is so since the JJA model considered here contains no dissipation mechanism 
which would arise e.g. in the presence of disorder or from a coupling of the phase degrees 
of freedom to normal electrons (Ohmic damping). The results of the numerical calculation 
for the regular part of the conductivity are shown in Fig.§. 

At T = the singular part vanishes while the regular part can be evaluated explicitly 
close to the critical point with the result 
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Note, that at the T = transition, where the gap in the response function vanishes ((5 = 0), 
there is no universal dc conductivity as in the two-dimensional caseta. Universality emerges, 
however, in a different context. We will now present the scaling analysis satisfied by (j{oj) 
in the vicinity of the quantum phase transition Ej = EJ^^, where the temperature obeys 
< ksT << Ej. The behavior of the conductivity in this regime can be understood in 
terms of a universal scaling function that depends on a variable which measures the distance 
of the superconducting ground state from criticality. In the quantum critical region we write 
the spherical constraint (§) in a terms of a low-temperature expansion 
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where Ki{x) is the MacDonald functionll^ {i.e. second modified Bessel function) and a 
Ec/ Ej. Solving Eq.(|l6]) for small parameter 5 we obtain: 






<rU-) = T-.l^^-^F (l^] ■ (17) 



The temperature ksT sets also the energy scale to measure the frequency via the ratio 
uj/kBT. It is therefore reasonable to introduce the dimensionless scaling variable X = u/kBT 
and finally explicitly write the F function as 




^;r ^ \3/2 



F{X) = -^X [ 1 - 4W — ^j coth(X), (18) 



where C2 = J^^ dxr]{x) / {3 — x)^/^. At criticality the power-law behavior of the model can 
be deduced from Eq.(^ by taking the momentum long wave-length limit. By generalizing 
the analysis to (i-spatial dimensions we find that the dynamical critical exponent z = 1, 
and that the correlation length exponent z/ = l/{d — 1), for d < du = 3, about the T = 
quantum critical point. At finite-temperatures, close to the quantum phase transition we 
obtain z = 1 and v = l/{d — 2), below the upper critical dimension du = 4, respectively (the 
anomalous dimension is in turn r/ = for all values of d). 

In conclusion, we have studied a 3-D quantum Josephson junction array model in the 
non-perturbative quantum spherical model approximation. We have explicitly calculated 
the phase diagram at zero temperature as well as the conductivity and its scaling properties 
a low temperatures about the quantum critical point. There are several problems left to 
consider in the future, like the role of disorder, dissipation, applied magnetic fields and the 
impact of anisotropy (relevant for layered high-temperature superconductors). 
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FIGURES 

FIG. 1. Charging energy parameter for the 3-D JJA as a function of the ratio of mutual Cm 
and self-capacitance Cg- 

FIG. 2. Phase diagram for a 3-D JJA in the parameter space defined by temperature T charging 
energy Eqc and the ratio of mutual- to self capacitance Cm/Cs of a single junction. The system is 
phase coherent in the region below the surface. 

FIG. 3. Real part of the dynamical conductivity at T = 0, for several values of the dimensionless 
gap parameter 6, that measures the distance from the critical point: 6 = 2, 6 = 0.5, 6 = 0.1 and 
6 = 0.05 (from the left to right). The arrow indicates the position of the singular Drude part 
(vanishing for T = 0). 
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